Quantum Confinement in Si and Ge Nanostructures 
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We apply perturbative effective mass theory as a broadly applicable theoretical model for quantum 
confinement (QC) in all Si and Ge nanostructures including quantum wells (QWs), wires (Q-wires) 
and dots (QDs). Within the limits of strong, medium, and weak QC, valence and conduction 
band edge energy levels (VBM and CBM) were calculated as a function of QD diameters, QW 
thicknesses and Q-wire diameters. Crystalline and amorphous quantum systems were considered 
separately. Calculated band edge levels with strong, medium and weak QC models were compared 
with experimental VBM and CBM reported from X-ray photoemission spectroscopy (XPS), X-ray 
absorption spectroscopy (XAS) or photoluminescence (PL). Experimentally, the dimensions of the 
nanostructures were determined directly, by transmission electron microscopy (TEM), or indirectly, 
by x-ray diffraction (XRD) or by XPS. We found that crystalline materials are best described by 
a medium confinement model, while amorphous materials exhibit strong confinement regardless 
of the dimensionality of the system. Our results indicate that spatial derealization of the hole in 
amorphous versus crystalline nanostructures is the important parameter determining the magnitude 
of the band gap expansion, or the strength of the quantum confinement. In addition, the effective 
masses of the electron and hole are discussed as a function of crystallinity and spatial confinement. 

PACS numbers: 73.21.-b,73.22.-f,78.67.-n,61.46.-w,81.07.-b 



I. INTRODUCTION 

Semiconductor nanostructures (NSs) exhibit increased 
oscillator strength due to electron hole wave function 
overlap, and band gap engineering due to the effect of 
quantum confinement (QC). Thus, materials like Si are a 
viable option for opto-electronics, photonics, and quan- 
tum computing. [1-3] QC is denned as the modification 
in the free particle dispersion relation as a function of 
a system's spatial dimension. [4] If a free electron is con- 
fined within a potential barrier, a shift in the band gap 
energy is observed, which is inversely proportional to the 
system size squared, in the effective mass approximation. 
As a result, the emitted photon energy is directly propor- 
tional to the gap energy (Eg). QC often manifests itself 
in optical experiments when the dimension of the sys- 
tem is systematically reduced and an increase in the ab- 
sorbed/emitted photon energy is measured correspond- 
ing to electron transitional states, i.e. in semiconductor 
NSs. 

For practical applications, utilizing QC effects in NSs 
requires an understanding of the band structure of a low- 
dimensional material, how the method of preparation ef- 
fects the final properties of the NS, and the kinetics/ 
dynamics of the absorption/emission process. The con- 
finement potential is determined by the alignment of the 
respective Fermi levels when a material of a Egi is sur- 
rounded by a material of a Eg2, with Egi < ^G2-[5] The 
preparation technique can introduce stress in the system, 
which changes the band gap energy. [6] For indirect gap 
materials phonon processes can effect the recombination 
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mechanism. [7] The lifetime associated with the recombi- 
nation event can be altered by the excitation power. [8] 
(For a review of general properties of low-dimensional 
structures, see Refs. 2, 4, and 9. For a discussion of 
other higher order effects in NSs, see Ref. 9.) This ar- 
ticle is concerned with the electron/hole recombination 
process in amorphous (a) versus crystalline (c) NSs with 
different dimensions. 

Several theoretical models (e.g. see Refs. 10-12) have 
been applied to NS; all models are empirical and no one 
model can model all semiconductor NSs. Since the pa- 
rameters of a NS system are dependent upon the prepa- 
ration method for a particular material, a comprehensive 
theoretical understanding must test along this dimension 
as well. In this article, we consider a relatively simple 
model of direct e-h recombination using a 'particle in a 
box' type model as a perturbation to the effective mass 
theory. We use no adjustable parameters [? ] and in- 
clude corrections to the model dependent on the prepara- 
tion method as known experimentally and/or computa- 
tionally when needed, thus achieving transparency in the 
physics involved. The only parameter tested in this work 
is the crystallinity, which is shown to effect the strength 
of confinement (defined in Sec. II), because of the differ- 
ent symmetry properties of the electron and hole. 

The model is applied to experimental results on crys- 
talline and amorphous Si and Ge NSs, including quan- 
tum wells (QWs), wires (Q-wires) and dots (QDs). Sys- 
tems of regular shape are chosen to ensure crystallinity 
is the primary parameter. For example, data obtained 
by van Buuren et al. [13] for high quality 'star-shaped' 
samples are difficult to analyse theoretically. Parameters 
relevant to a particular system are discussed and energy 
corrections are given when needed. Briefly, we compare a 



few theoretical models with experiment, thus, illustrat- 
ing the need to categorically understand experimental 
parameters. Results are discussed and a mechanism for 
the differences between the strength of confinement in 
the amorphous and crystalline system is proposed. 



II. THEORY 

In this work, we use the effective mass approximation 
(EMA) based on the Bloch periodic function. The essen- 
tial features of the model are discussed below. 

The Bohr radius of an electron (e) , hole (h) or exciton 
(X) is given by, in SI units: 



a e (h)(x) 
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is the effective mass of the e, h or X, respec- 
tively, e is the electric charge and e is the dielectric con- 
stant. Depending on the e or h effective mass, the X-Bohr 
radius is 4.5 nm for Si and 24 nm for Gc. The Bohr radius 
defines the spatial dimension of the particles, which de- 
termines the range of sizes for which QC can be observed. 
We define three regimes of confinement here: [4] 

• Weak confinement: When the dimension of the sys- 
tem is much larger than a e and a^. In this situa- 
tion, the appropriate mass in the kinetic term is 
M = to* + m* h . The energy term is dominated by 
the Coulomb energy. 

• Medium confinement: When the dimension of the 
system is much smaller than a e , but larger than 
ah, then only electrons will experience confinement. 
The relevant mass is simply to* for the kinetic term. 
Most materials belong to this class. 

• Strong confinement: When the dimension of the 
system is much smaller than a e and a^. Here both 
electrons and holes experience confinement and the 
relevant mass is the reduced mass, fx, with - = 
-=r + -V . In this regime, the Coulomb term is small 
and can generally be treated as a perturbation. 

Below we will use the terms 'weak,' 'medium' and 'strong' 
to refer to the different regimes of confinement discussed 
above. 

Si and Ge are both indirect gap materials, meaning 
that, in principle, phonon scattering events are essential 
to maintain momentum and energy conservation during 
a radiative event. This situation is true in the case of 
a bulk material; however, as the dimension of the sys- 
tem is reduced, the uncertainty in the momentum k vec- 
tor is increased. Therefore, it is possible to break the 
k selection rules making the Eq 'pseudo-direct,' allow- 
ing for direct e-h recombination. [14] The length scale at 
which this 'pseudo-direct' phenomenon becomes impor- 
tant is typically less than a few nanometres. [12, 15, 16] 



This length scale corresponds to the systems considered 
here; therefore, theoretically it is valid to assume direct 
e-h recombination without phonon-assistance. 

In the 'particle in a box' model the bulk Eq is taken 
as the ground state energy. The effect of reduced di- 
mension is considered as a perturbation to the bulk Eg- 
Therefore, we consider the general field Hamiltonian for 
a system of Coulombic interacting particles given by (de- 
tails are given in Rcf. 17): 



H = J d 3 r^(r) f^V 2 ) m + 



1 [ d 3 rd 3 f^(r)^(f) 6 ,M r)^{r), (1) 

2 J 47re|r — r\ 

where ip(r) is the field operator, to is the mass of the 
electron or hole, e is the dielectric constant of the sur- 
rounding medium and e is the electric charge. We do not 
consider the spin-orbit interaction here, because the fine 
structure is negligible at the energies considered here. 

The field operators are expanded in a two-band model 
for the conduction band C and the valence band V as: 



^( r ) = ^2 a k,i ( Pks( r ) (* ^ C, V), 



(2) 



where k represents a summation over momentum states. 
The ifik.i(r) basis set in Eq. (2) is expanded to reflect the 
use of an infinite confinement potential with a Bloch basis 
Uk,i- Infinite confinement is a reasonable assumption for 
the systems we are considering, because the matrix mate- 
rial has a Eq several eV higher than the nano-structure; 
however, we can not discuss hopping or other such higher 
order effects. Bloch states reflect the periodic nature of 
the crystal (Luttinger-Kohn representation), while the 
boundary conditions of a NS do not reflect this same pe- 
riodicity. However, in many NSs the transitions we are 
interested in happen near the Brillouin zone centre, e.g. 
the T-point. This statement may not be strictly true in 
the case of weak confinement, because fc-selection rules 
are not as strongly broken as in the case of strong confine- 
ment. Nonetheless, kp perturbation theory considers ex- 
pansions about the Brillouin zone minimum, k D . There- 
fore, we may justify the use of Bloch states through the 
use of the slowly varying wave approximation whereby 
only the k o =0 states are retained. 

For indirect gap materials the exciton is Wannier-likc, 
in the limit k <C — (a c is the lattice spacing) and we 
can drop the exchange term, which goes to zero quickly. 
Equation (1) is solved in the exciton basis using the state 
$ defined as an e-h pair above the ground state, $o, as: 

$ = 'LkikzCkikAX^V, & $V = h 3 b ki ---b kN $ , 

where a k (bk) refers to electrons (holes) in the conduction 
(valence) band. Expanding in low lying fc-states near the 
band edge, we solve Eq{D) = ($| H which gives the 
variation of gap energy with nano-structure size. 

For the mass terms in Eq. (1), we use the effective 
masses calculated using the density of states. [18] The ef- 
fective mass is related to the parabolicity of the band 



TABLE I. Parameter A given in Eq. (3) for 3D, 2D, ID 
confinement and for AEcbm, AEvbm- 
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structure, which is not expected to change in a nano- 
structure compared to a bulk material at the T-point. 
Therefore, we assume the effective mass from the bulk 
system. For Si the effective masses at room temperature 
are: m c — > m* = f.08m o and my — > m v = 0.57m o . 
For Ge the effective masses arc: m c — > m* = 0.56m o 
and my —5- m v = 0.29m o . These definitions yield the 
equation: 

E Gap {D) =E Gap {<x) + -^eV-mn 2 . (3) 

Ego, p (oo) is the band gap of the bulk material and D 
represents the QD diameter, the QW thickness or the 
Q-Wire diameter in what follows. The calculation was 
carried out for confinement in ID, 2D with cylindrical 
coordinates and 3D with spherical coordinates. The pa- 
rameter A is given for Si and Ge in the strong, medium 
and weak confinement regimes in Table I. The change in 
energy of the CBM (AEcbm) due to QC is labelled as 
'medium confinement' in Table I, because a AEcbm is 
equivalent to QC of the electron only as defined by our 
model, where only the electron mass is considered in Eq. 
(1). The change in energy of the VBM (AEvbm) due to 
QC is also listed in Table I, which is calculated by con- 
sidering confinement of the hole only, where only the hole 
mass is considered in Eq. (1). The other fixed parameter 
is the appropriate Eg (00) of the bulk system and one 
could argue for the use of a renormalized effective mass 
with dimension of the system, which is discussed in Sec. 
V. 

Finally, it is important to note that theoretical mod- 
elling can be further complicated by the accuracy of NS 
size determination. Transmission electron microscopy 
(TEM) is the direct method to determine NS size; how- 
ever, if the contrast between the matrix and the nano- 
structure is poor, then the size uncertainty can be on 
the order of 1 nm.[19] Indirect size determinations can 
be used as well, such as with x-ray diffraction (XRD) 
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FIG. 1. Disordered Si-QW data and theoretical fit. Experi- 
mental data from Ref. 27. Theoretical fit using ^4=0.89 and 
Ecap(oo) = 1.6 eV in Eq. (3). NB: The CBM shift is offset 
by the E G a P (oo). 



[20] or x-ray photo-emission spectroscopy (XPS).[21] Fur- 
thermore, QC in Ge has been a greater challenge for 
researchers to observe than in Si, because of the ten- 
dency to form defects, interfacial mixing and sub-oxide 
states. [22-25] Therefore, only limited results on Ge are 
discussed here. However, there is recent progress in this 
area, showing very promising results. [26] 



III. EXPERIMENT 

We cite the results of several experimental works in- 
cluding our own from the University of Western Ontario 
and from the National Research Council Ottawa, in Sec. 
IV. The essential features of each experiment are given 
here. The details of the experiments can be found in the 
references provided. 



IV. RESULTS 
A. Silicon 

1. Quantum Well 

Si/Si02 superlattice Si-QWs have been grown using 
molecular beam epitaxy, determined to be disordered 
via Raman scattering measurements, and their thickness 
found using TEM and XRD. [27, 28] The change in the va- 
lence band maximum (VBM) and conduction band min- 
imum (CBM) position was measured using XPS and Si 
L2,3 edge absorption spectroscopy, respectively, and room 
temperature PL spectroscopy was measured. Fig. 1 plots 
the model predictions with the experimental data. 

In Ref. 27 the authors used a fitting procedure accord- 
ing to the effective mass theory for the AE V bm(gbm)i 
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FIG. 2. Crystalline Si-QW data and theoretical fit. Experi- 
mental data from Ref. 21. Experimental PL data from Ref. 
29. Theoretical fit using A=0.35 and E Gap (oo) = 1.12 eV in 
Eq. (3). NB: The CBM shift is offset by the E Gap (oo). 



resulting in AE VB m = -0.5/D 2 and AE C bm = 0.7/£> 2 , 
where D is the thickness of the QW. Our model predicts 
AEvbm = -0M/D 2 and AE CBM = 0.35/Z? 2 . The 
trend for AEcbm is more accurately given in Ref. 27. 
In Ref. 28, the change in Eq was fitted with A = 0.7 
and EGap(oo)—1.6 eV, as in Eq. (3). The fit also de- 
termined the effective mass to be m£( e ) ~ 1. The model 
uses EGap(oo)=1.6 eV to fit the experimental PL data 
well when employing the curve for strong confinement 
with A = 0.89. 

Next we look at c-Si/Si02 QWs fabricated by chemical 
and thermal processing of silicon-on-insulator wafers. [21] 
The same methods described above were used to deter- 
mine experimentally the AE VBM (gbm) and the change 
in the gap energy including the total electron yield for 
a better signal to noise ratio. The thickness of the Si 
layer was determined by XPS using a mean free path 
in Si of ~1.6 nm. Note that a thickness of 0.5 nm corre- 
sponds to a single unit cell of Si. Therefore, experimental 
data below w 1 nm should be treated with caution. In 
a parallel study, these c-Si/Si02 QWs were investigated 
optically. [29] 

Fig. 2 compares experimental measurements and the 
model results for c-Si-QWs. The Eg(oo) in the model is 
1.12 eV and the AEvbm is not significant below 1.5 nm. 
The AEcsM) AEgbm+vbm , and the experimental PL 
are all well fitted by the curve for medium confinement, 
with A = 0.35. In Ref. 29 it was found that there is a sec- 
ond PL peak fixed with respect to the Si layer thickness 
at 1.8 eV. This second peak was associated with inter- 
face states. Therefore, we can assign the experimental PL 
data in Fig. 2 with direct e-h recombination modelled by 
medium confinement. 



FIG. 3. Crystalline and amorphous Si-QD data and theo- 
retical fit. 'Expt. Ion-Implantion Si02' refers to crystalline 
Si QDs embedded in Si02 from Ref. 30. 'Expt. microwave 
plasma decomposition (MPD) Si02' refers to crystalline Si 
QDs embedded in Si02 from Ref. 20. 'Expt. plasma en- 
hanced chemical vapour deposition (PCVD) SiN' refers to 
amorphous Si QDs embedded in SiN from Ref. 34. Theoret- 
ical fit using A=3.57 and 1.39 and EG ap (oo) = 1.12 or 1.56 
eV (as labeled) in Eq. (3). NB: The absorption data is offset 
by the E Gap (oo). 



2. QDs 

First we consider Si QDs formed by ion implantation 
in Si02 films, followed by high-temperature annealing 
in N2 and forming gas. [30] Ref. 30 reports the QD di- 
ameter and crystalline structure observed by TEM, and 
room temperature PL measurements. TEM data show a 
Gaussian distribution in the Si-QD diameter with depth, 
resulting in a stretched exponential PL dynamic. [30, 31] 

We compare ion-implanted Si-QDs with Si QDs in a 
Si02 matrix prepared by microwave plasma decompo- 
sition (MPD) creating ultrafine and densely packed Si 
QDs [20] (implying that tunnelling effects are important 
here [32]). The crystallinity and size was determined by 
TEM imaging and XRD, respectively. In Ref. 20, the au- 
thors note that PL was not observed unless the Si QDs 
were oxidized, implying that surface bonds were passi- 
vated with suboxide states eventually forming a surround 
Si02 matrix. 

Fig. 3 shows the experimental PL data for ion- 
implantated and MPD Si QDs together with our cal- 
culated curves for strong and medium confinement. 
Above 3 nm both sets of experimental data follow 
closely the model of strong confinement with A=3.57 and 
Eg(oo)=1A2 eV. This indicates that for sample diame- 
ters larger than this size tunnelling effects are significant, 
implying a de-localization of carrier states. Iacona et al. 
measured a similar trend for experimental PL data. [33] 
Below 3 nm, when QC effects are particularly strong, 
the ion-implantation data follows the curve for medium 
confinement, with yl=1.39. 



Next we consider a-Si QDs embedded in a SiN 
matrix. [34] The Si QDs were fabricated using plasma en- 
hanced chemical vapour deposition. The size and amor- 
phous structure were measured using TEM and the PL 
was taken at room temperature. Absorption data was 
taken by ultraviolet-visible absorption spectroscopy. The 
value for the bulk band gap given by the authors is 1.56 
eV, which is obtained via a fitting procedure. This value 
is known to vary between 1.5— >1. 6 cV, for Si samples 
prepared similarly. [34] 

We can see in Fig. 3 that the experimental data 
for absorption and PL of a-Si QDs embedded in SiN 
lies between the curve for medium (A=1.39) and strong 
(A=3.57) confinement, with £' C ; Q p(oo)=1.56. Using a fit- 
ting procedure, the authors of Rcf. 34 found A=2.40. 
The authors further conclude that by observing the fact 
that the experimental absorption data lies close to the 
PL data, one can conclude that the PL data for these 
samples is a good measure of the actual change in the 
Eg(D).[3A] Notice that this situation is similar to that 
observed for Si-QWs (see Fig. 1 and 2). 



3. Quantum Wires 

Due to inherent complications in the fabrication pro- 
cess of Si or Ge wires with a diameter below the Bohr 
radius, few studies on QC in nano- wires exist and we are 
only able to report on c-Si-Q-wires. On the other hand, 
por-Si studies are widely cited in the literature. With 
suitable control of the etchant, por-Si QDs can become 
elongated, [35] thus breaking confinement in one direction 
implying they are more wire-like; a detailed discussion is 
provided in Ref. 36. In this case, they are called pseudo- 
por-Si-QDs or in the case they behave like interconnected 
dots, spherites.[37] 

Anodically grown por-Si samples were prepared by 
Schuppler et al.[38] X-ray absorption measurements de- 
termined the structures to be closer to c-Si than to a-Si. 
TEM was used to determine the size and PL measure- 
ments were performed at room temperature. The por-Si 
structures are said to be H-passivated and O-free; how- 
ever, samples were exposed to air. 

Si Q-wires were produced by Ma et al. using an oxide- 
assisted growth method with SiO powders. [39] Subse- 
quently, the wires were cleaned with HF to remove the 
oxide, thus forming a H-terminated surface. Scanning 
tunnelling microscopy was used to determine the diame- 
ter of the wires. The formation of SiH2 and SiH3 was ob- 
served on the facets of the Q-wires, which was attributed 
to bending stresses in the wires. The energy gap was de- 
termined using scanning tunnelling spectroscopy, which 
also indicated doping levels in the wires as seen by an 
asymmetrical shift of the Eq around V. 

The experimental data from Ma et al. and Schuppler 
et al. can be seen in Fig. 4. Below 3 nm the experimental 
data from Schuppler et al. ('por-Si Wire PL') lie close 
to the curve for 2D strong confinement with A— 2.09 and 
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FIG. 4. Crystalline Si Q-wire and QD data and theoreti- 
cal fit. Experimental por-Si wire data from Ref. 38. Ex- 
perimental Si Wire data from Ref. 39, using scanning tun- 
nelling spectroscopy (STS). Experimental por-Si QD data 
from Ref. 37 and 41. Theoretical fit using A=1.39 and 0.81 
and EGap(oo) = 1.12 eV in Eq. (3). NB: The absorption data 
is offset by the EG ap (oo). NB: 'QD' here refers to spheroids. 

EG a p(oo)=1.12 eV. Notice that the experimental data 
also lie close to the curve for 3D medium confinement 
with A=1.39. This observation may be a reflection of the 
idea that these structures are between dots and wires. 
On the other hand, the data from Ma et al. lie close 
to the curve for 3D strong confinement, using the same 
EGap(oo) and A=3.57. We also note that recently Si-Q- 
wires have been produced [40] with results nearly identical 
to those of Ma et al. 

Experimental data on pseudo-por-Si-QDs for both ab- 
sorption and PL are taken from Ref. 37 and 41. Raman 
and TEM measurements were used to determine the size 
and the 'sphcritc' nature of the samples, respectively. PL 
measurements were performed at room temperature and 
at 4.2K, with very little difference in the two measure- 
ments. Optical absorption was performed at room tem- 
perature. It is also noted in Ref. 41 that, for por-Si, 
interface states and phonon events are significant. Fig. 4 
shows the PL and absorption experimental data for por- 
Si-QDs. Here the experimental data are modelled by the 
curve for 3D strong confinement, with ^4=3. 57 and the 
same gap energy as above. Compared to absorption and 
PL data for a-Si-QDs in Fig. 3 and the Si-QWs in Fig. 
1 and 2, there is a significant shift between the absorp- 
tion data and the PL data, indicating a Stokes shift in 
the emission. [41] Furthermore, as noted in Ref. 37, the 
experimental PL data are nearly identical to Takagi et 
al., shown in Fig. 3. 

B. Germanium 

The first observation of QC in Ge was by Takeok et 
al.[42] In this study, they produced Ge QDs using an 
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FIG. 5. Crystalline Ge-QD data and theoretical fit. Exper- 
imental Ge QDs data from Ref. 42. Experimental Ge ab- 
sorption data from Ref. 43. Theoretical fit using A=2.69 and 
E Gap (oo) = 0.66 eV in Eq. (3). 



rf co-sputtering method followed by thermal annealing. 
The size of the Ge QDs was controlled by varying the ini- 
tial Ge concentration and was later determined by TEM 
imaging, which also showed that the Ge QDs were highly 
crystalline. PL was performed at room temperature. 

In a more recent study, Ge QDs were produced by con- 
densation out of the gas phase onto a Si substrate cleaned 
by HF.[43] The Ge QDs were determined to be in the bulk 
diamond crystalline phase. X-ray absorption (XAS) data 
were taken and can be seen in Fig. 5. XAS excites the Ge 
2p electron into the conduction band; therefore, the re- 
searchers obtained data for the change in the conduction 
band. 

The experimental data from Refs. 42 and 43 are pre- 
sented in Fig. 5. Note that the absorption data are 
obtained by shifting the Bostedt et. al. data by the 
Ec(oo) of Ge at 0.66 cV. Further note that above 3 nm 
there is a nearly identical departure from the medium 
confinement curve into strong confinement as was seen 
with Si-QDs in Fig. 3. In general, both sets of experi- 
mental data are well modelled by the curve for medium 
confinement with A=2.69 and Ec a p(oo) = 0.66 eV. For 
the smaller sizes (below 2.5 nm) the behaviour appears 
to deviate from medium confinement. This result may 
be because for the smaller sizes the authors only esti- 
mated the sizes. [42] In Ref. 43 the Ge-QD diameter was 
determined using atomic force microscopy, which can po- 
tentially give a larger uncertainty in determining the size 
of the dot. [44] Therefore, if the QDs are not symmetric 
then the diameter measurements could be inaccurate. 



V. DISCUSSION 

We start our analysis by giving a justification of the 
EMA, while highlighting some of the limitations. Exten- 



sive arguments appear in the literature concerning the va- 
lidity of the EMA and its k • p generalization. On the one 
hand, it is argued and demonstrated that the EMA over- 
estimates the -Eg [4, 10]; however, Sec. IV demonstrated 
that in some cases the EMA can underestimate the Eq- 
In part, this is because due to QC the parabolic nature of 
the bands is possibly removed. Another complication can 
arise from the fact that the envelope functions may not 
be slowly-varying over the unit cell, which is essentially 
complicated by the boundary conditions. A central prob- 
lem for EMA is in its applicability to a-materials, because 
it is based on the assumption of translational symmetry. 
Street has argued that while it is strictly not justified in 
the a-system, due to nonspccifically-dcfincd k vectors, it 
is still widely used albeit with differing assumptions. [45] 
We will discuss further the application of the EMA to 
both amorphous and nanostructured-systems below. 

On the other hand, it has been argued by See et al. 
that the EMA is well justified and produces agreement 
with the tight-binding method. [46] Such arguments re- 
side in the fact that it is not clear what all the relevant 
parameters are in a nano-structured system of a particu- 
lar material. In general, the boundary conditions of the 
system become very important, which is a problem for 
all theories. [47] If the Fourier components of the enve- 
lope function are centred around the the Brillouin zone 
centre, then envelope functions can be justified. In addi- 
tion, this justification has been extended to consider that 
if the interface is defect free then the EMA is justified. [47] 
Other considered corrections to the EMA use a fourth or- 
der term in k.[4] The advantage of the EMA is that it 
is straightforward in its application, thus allowing one 
to highlight key features of individual systems. Pertur- 
bations in the NS system arc naturally treated in the 
kp method and defect states easily calculated[18]. Com- 
pared to empirical methods, [10, 48] which produce a di- 
mensional dependence of D -139 , the EMA has the units 
D~ 2 (see (3)). In addition, it has been shown that the 
k • p Hamiltonian can be made to reproduce multiband 
coupling effects and the correct symmetry of the QD.[49] 

To emphasize the importance of accurately parametriz- 
ing the preparation method, we compare our results with 
a few theoretical models with respect to experiment. In 
the work by Bulutay[50], the variation in the F,q(D) is 
calculated using an atomistic pscudopotential method. 
The result is given in Fig. 4 of Ref. 50 and is repro- 
duced here in Fig. 6. References for the works listed in 
figure caption are given in Ref. 50. The top curve in 
Fig. 6 is for Si, which we compare with our Fig. 3 here 
and the bottom curve for Ge, which we compare with 
our Fig. 5 here. It is clear that the band gaps shown in 
Fig. 6 are consistently larger than what we present in 
this manuscript. This result is easy to explain. 

In the case of Si, the experimental data in Fig. 6 is 
from Furukawa et al. In this work, they produce Si QDs 
using magnetron rf sputtering. It is demonstrated that 
the QDs are surrounding by H and composed of Si:H. 
The incorporation of H in the QDs causes an increase 
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FIG. 6. Comparison of the Eg as a function of QD diameter 
for Si and Ge. Reprinted figure with permission from Ref. 50. 
Copyright (2007) by the American Physical Society. Refer- 
ences for the works listed in figure caption are given in Ref. 
50. 



in the Eg (see Fig. 1 of Ref. [51]). Furthermore, the 
Raman peak of the QDs is measured to be 514 cm -1 , as 
opposed to the bulk value at 520 cm" 1 , which indicates 
that the system is under stress. To verify this claim, the 
authors measure a 2% extension in the bond length us- 
ing x-ray diffraction. Thus, the QDs are under tensile 
stress, which increases the Eg [52], see below for more 
detail. Finally, Furukawa et al. explain that the origin of 
the stress is through the incorporation of the above men- 
tioned H in the Si QD lattice through a plasma-assisted 
crystallization process. Whereas, if H was acting only to 
passivate the dangling bonds there would be no change 
in the Eg. [53] Therefore, the experimental results of Fu- 
rukawa et al. are higher than what is presented here (Fig. 
3) because of an additional stress component, which in- 
crease that Eg beyond that of QC alone and which is 
introduced because of the preparation method. 

The remaining results in Fig. 6 for Si are theoretical 
results. The work by Ogiit et al. uses a real-space pseu- 
dopotential method. Vasiliev et al. uses linear-response 
within the time-dependent local density approximation. 
Garoufalis uses time-dependent density functional the- 
ory. Including the work of Bulutay, all four methods give 
approximately the same result. However, looking at Bu- 



lutay, there is no explicit inclusion of a stress component 
in the Hamiltonian, instead it is implicitly fitted in the 
pseudo-potential, while the other three methods ignore 
stress altogether. Therefore, these four methods do not 
explicitly consider the experimental details, instead they 
are fitted to experiment. The empirical nature of the 
theoretical methods can be further seen when comparing 
with similar pscudopotcntial calculations that produce 
different results from those shown here. [10] 

Considering the results for Ge in Fig. 6, the situation 
is essentially the same as for Si above. The experimental 
data from Kanemitsu et al. is associated with defect PL 
only, making it beyond the scope of this paper. Niquet et 
al. uses an sp 3 tight binding method, while Tsolakidis et 
al. uses time-dependent density functional theory in the 
adiabatic local density approximation. The experimental 
data from Takeoka et al. is fitted in this manuscript (Fig. 
5) and not in the references of Fig. 6. Furthermore, these 
calculations for Ge are similar to Si, which implies the 
material properties are not being properly accounted for 
theoretically. Thus, it is necessary to quantify each term 
in the Hamiltonian according to the preparation method. 

Finally, we comment on the relevant energy scales for 
the experiments considered above. The electron and 
hole can form a hydrogenic or positronium-like exciton, 
a bound state of the constituent particles, thus modify- 
ing the photon energy during the recombination event by 
the Coulomb interaction between the electron and hole. 
The Coulomb energy is on the order of hundreds of meV 
(— > l/R, i?=NS dimension), the exchange energy is on 
the order of 0.1 meV (— >• 1/i? 3 ), while the gap energy 
is on the order of several eV. Due to the large number 
of competing parameters in any real system, the exact 
value of the above parameters is not known, and these 
are important for precision control of a device. 

To summarize the comparisons made in Sec. IV, we 
first consider the relationship between experimental ab- 
sorption and PL data. In the case of disordercd-Si-QWs 
(Fig. 1), c-Si-QWs (Fig. 2) and a-Si-QDs in SiN (Fig. 3) 
the absorption curve follows closely with the PL. As men- 
tioned in Sec. IV A 2, this result indicates that the PL 
measurement is an accurate measure of Ec{D). Further- 
more, in the case of Si-QWs the VBM does not change 
significantly. Therefore, we conclude that the model de- 
pendence between these three systems does not lie in the 
change in the VBM. 

Considering the absorption data from por-Si-QDs (Fig. 
4) , there is a significant shift between the absorption data 
and PL data, which was noted in Sec. IV A 3. In addi- 
tion, the por-Si QD data are nearly identical to the MPD 
Si-QDs (Fig. 3), which indicates that these systems are 
structurally similar with similar decay dynamics. In the 
case of por-Si it has been found that this system is under 
tensile stress. [54] Tensile stress, which is a function of the 
thickness of oxide, is known to increase the band gap. [52] 
It is known that the surrounding oxide has a strong effect 
on the resulting PL in por-Si. [55] The resulting Si-O-Si 
bonds due to the oxidation process place large stresses 



on the por-Si crystallites. In addition, it has been shown 
that the dominant PL comes from surface states. [41] At 
the surface or interface states, it has been shown that 
band bending on the order of 0.2— >0.3 eV can occur. [56] 
Such a shift in energy corresponds with the discrepancy 
shown in Fig. 3 and 4. 

For the c-Si-Q-wires measured by STS (Fig. 4), the 
data are modelled by strong confinement. This is be- 
cause of the stresses observed in the system and possibly 
because of the doping; both are factors that do change the 
nature of electronic structure. In Sec. IV A 3, we men- 
tioned that these structures experience bending stress, 
which has a tensile component. Furthermore, Fig. 4 il- 
lustrates that c-Si-Q- wires are identical in energy to por- 
Si; therefore, the analysis of these systems is similar. By 
contrast, the por-Si wire PL data (Fig. 4) behaves more 
wire-like, which may be a result of the fact the authors 
took care to minimize oxygen exposure (see Sec. IV A 3). 

From Fig. 3 and 5, both ion implanted Si-QDs and Ge- 
QDs have the same behaviour above 3 nm. They lie close 
to the curve for strong confinement, similar to the case 
of por-Si, indicating that possible stresses or interface 
states are important in this regime. Ge is known to ex- 
perience stress in a Si02 matrix. [57] Tensile stress can be 
relieved depending on the nature of the interface bonds 
and the surface to volume ratio of Si:Si02-[52] In the 
work of Rcf. 17 it was found from Raman spectroscopy 
that ion-implanted QDs arc not under stress for diame- 
ters smaller than 3 nm. Therefore, c-Si-QDs produced 
by ion implantation and c-Ge-QDs are well modelled by 
medium confinement below 3 nm. 

Finally, a-Si-QDs in SiN (Fig. 3) lie between medium 
and strong confinement (see Sec. IV A 2) . SiN has a band 
gap of 5.3 eV versus Si02 at 9.2 eV, which allows for 
tunnelling of carrier states. [34] More importantly, if we 
consider the nucleation process during thermal anneal- 
ing and consider the bond enthalpies for diatomic species 
(SiN at 470 kJ/mol and SiO at 799 kJ/mol), it is easier 
to break SiN bonds, thus allowing for a greater degree 
of intermixing at the QD-matrix interface. Therefore, a 
SiN matrix acts more like a finite potential barrier, which 
lowers the gap energy from the infinite case. A numeri- 
cal computation indicates that the difference between the 
case of finite versus infinite confinement potential is be- 
tween 10% and 15% depending on the size of the system. 
This difference exactly corresponds with the difference 
we see in Fig. 3. Therefore, we conclude that a-Si-QDs 
in SiN are well modelled by strong confinement. 

From the results above and considering modifications 
that must be made to our model to account for non-direct 
e-h recombination phenomena, it is clear that strong con- 
finement describes a-materials and medium confinement 
describes c-materials. Therefore, since QC of a particle 
is a function of the derealization of that particle with 
respect to the dimension of the system, we need to ac- 
count for the fact that the hole becomes more delocalized 
in the a-system than in the c-system. This fact may or 
may not be seen as a shift in the VBM. As noted above, 



disordered- Si- QWs, c-Si-QWs and a-Si-QDs in SiN all do 
not show a large variation in the VBM. 

A mechanism for pinning of the hole states in c-Si-QDs 
was discussed in the work of Sa'ar et al. as a function 
of the hole coupling with vibrons.[58] However, this phe- 
nomenon does not account for the fact that the hole be- 
comes more delocalized in the a-system, it is well known 
that band-tail states play a very important role in the 
band structure of a-materials, even though the popula- 
tion density is relatively low. [45] Kanemitsu et al. (and 
Refs. within), report the experimental observation that 
the band-tail states become strongly delocalized in the 
a-system, while the hole remains relatively localized in 
the c-systcm.[59] This observation accounts for what is 
observed in this work. 

Another critical factor to discuss is the effective mass 
concept, particularly in the a-systcm. Recall from Sec. 
II, the bulk effective mass is used in the calculations. 
It is possible that this parameter is not well-justified in 
the a-system[45] and is simply not valid in the nano- 
structured system, in the worst possible case, or it is 
size-dependent. [60-62] 

The electron (e) and hole (h) interact differently with 
the atomic structure. The s-like electron has rj? symme- 
try; the p-like hole is contained in r®. Therefore, holes 
interact more strongly with the acoustical lattice vibra- 
tions. The electron has approximately twice the mass of 
the hole, which is dependent on the gap energy. Hence 
the crystallinity will effect the properties of the particles 
differently and recombination events are dependent on 
such properties. 

The a-system has typically 80% [63] of the density of 
the c-system, while disordered Si generally refers to a 
density ss 98% [64] of that of c-Si, and these values can 
vary widely based on the preparation method. [65] There- 
fore, short and medium range structural order does re- 
main in both of these systems. Although the long-range 
order is not well-defined in the a-system along with the k- 
vectors, alternative approaches to this concept have been 
extensively presented. In an earlier work, Kivclson et al. 
defined an alternative approach to this concept. [66] They 
formulated the assumptions (i) the structure of the solid 
can be approximated by a rigid continuous random net- 
work that is homogeneous on the scale of the slowly vary- 
ing envelope, and (ii) the band can be measured by a set 
of linearly independent orbits, which are not necessarily 
orthogonal. Furthermore, Kivclson et al. used a tight- 
binding approach with approximate eigenvalues to obtain 
the effective mass Hamiltonian. In another approach[67], 
Singh looked at the effective masses in the extended and 
tail-states around the mobility edge directly using a real- 
space formulation. The electron energy eigenvalues are 
given in terms of the probability amplitude, which can- 
not be defined as in the case of a c-material in terms of 
k- vectors. Instead, the probability amplitude is defined 
as [67]: 



Cn = N~ 1/2 cxp(is e ■ 1), with s e (E) = 



where Ec defines the mobility edge; therefore, the ef- 
fective mass is defined above the mobility edge in the 
extended states and is imaginary in the tail states. In 
either approach described here, the result is that the 
effective mass calculated is lower than in the bulk sys- 
tem. This observation implies that the Bohr radius of 
the hole in the a-material is larger than in the c-material 
and hence the hole is more delocalized in the a-material, 
thus the observed strong confinement. It is clear that 
this is the dominant mechanism for strong confinement 
in the amorphous system, since the pinning discussion 
above does not describe all the systems considered here. 
The relative magnitude of the two mechanisms needs fur- 
ther analysis. 

The size dependence of the effective mass in c-systems 
is reported in Refs. 60-62. Experimentally, the effective 
mass is reported to decrease with size in Ref. 61 and 
62. In one theoretical report, the hole effective mass in- 
creases, while the electron effective mass decreases. [60] 
The magnitude of change in the effective mass is roughly 
the same for the electron and the hole, and consider- 
ing the effective mass of the electron in the bulk sys- 
tem is roughly twice that of the hole, it is not likely 
that the change will be within experimental resolution. 
Overall, the effective mass in the a-system and in the 
nano-structured system is understood to decrease, but 
the magnitude of the decrease is unclear. Therefore, in 
terms of the calculations presented here, if the effective 
mass is lowered than we should expect to see an increase 
in the calculated Eq and hence our curves will shift up- 
wards. However, we would also expect to see an increase 
in Eq from the experimental results. Since the exact 
magnitude of these changes is not known it is difficult 
to evaluate the error incurred by using the bulk effective 
mass. 

This issue of the correct effective mass is more poignant 
when considering the AEc BM (VBj\/) . In this work, 
AEysAf > AEcbm, which is understood, because the 
effective mass of the hole is smaller than the electron. 
However, experiment consistently shows the opposite ef- 
fect, see Fig. 1,2 and see Ref. 68 and 69. This ob- 
servation implies that experiment is measuring a larger 
decrease in the electron effective mass than the hole, or 
possibly a relative increase in the hole mass compared to 
the electron. This observation is nearly consistent with 
Ref. 60, where they predict a nearly symmetric change. 
Furthermore, recall that experiment reports a decrease 
in the electron effective mass. [61, 62] 

Therefore, the decrease in the electron effective mass 
and increase in the hole effective mass is consistent for 
the crystalline system with our observation of medium 
confinement, because the hole is more spatially localized. 
In being consistent with experiment, we drop the hole 
contribution for the crystalline system in the ideal ap- 
proximation, because this term is not as significant as 



the electron according to the AEcbm(vbm) measure- 
ments, described above. Although, there may still be 
a slight hole contribution in this ideal approximation, 
which needs further study. In addition, in our theoretical 
modelling, we have consistent results for strong confine- 
ment in the amorphous samples, because both the elec- 
tron and hole effective mass decrease implying confine- 
ment of both, due to spatial de-localization. Although, 
the relative contribution from the electron versus the hole 
is not clear and needs further study. Furthermore, it is 
clear that the effective mass prediction for AE C 7bm(vbm) 
is not correct, unless a renormalized effective mass is used 
according to system dimension. These results are a clear 
indication that the use of the bulk effective mass is only a 
first order approximation. Nevertheless, very good over- 
all agreement is obtained between experiment and a the- 
ory with essentially no adjustable parameters for both Si 
and Ge nanostructures. 



VI. CONCLUSION 

We have studied the effect of confinement dimensions 
and crystallinity on the magnitude of the band gap ex- 
pansion (as a function of decreasing size) in group IV 
semiconductor NSs (quantum wells (QWs), wires (Q- 
wircs) and dots (QDs). Medium and strong confinement 
models provide the best fit to experimental results; more- 
over crystalline materials exhibit medium confinement, 
while amorphous materials exhibit strong confinement 
regardless of the confinement dimensions of the system. 
This difference in confinement strength was explained 
by considering the extent of spatial derealization of the 
hole. A possible explanation is hole pinning due to cou- 
pling with the vibronic states. [58] It has previously been 
reported [59] that band tail states become strongly delo- 
calized in the amorphous system compared to the crys- 
talline system. This hole derealization would partially 
account for the trends observed in our work. The con- 
cept of the effective mass was reviewed for the amor- 
phous system. We argue that the effective mass can still 
be defined in the amorphous material around the mobil- 
ity edge. [66, 67] A lower value of the effective mass is 
reported for the amorphous system, which accounts for 
the trends observed in our work, while the hole mass in- 
creases and the electron mass decreases as a function of 
spatial confinement. [60-62] With the diminished effective 
mass (the absolute value of this change is not possible to 
estimate, and more work is needed in this area), we ex- 
pect an increase in Ec ap , and our calculated curves of 
energy versus diameter will be shifted upwards. However 
the general trends observed in this work will remain the 
same. 
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